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Many anomalies of water have been discovered. Some concern the static thermodynamic properties 
of water, such as the increase of compressibility and specific heat upon decreasing the temperature. 
Others concern the dynamic properties, such as the breakdown of the Stokes-Einstein relation in 
supercooled water ^ and the non-Arrhenius to Arrhenius dynamic crossover 

The thermal conductivity A is the proportionality constant between heat flux and temperature 
gradient that causes the flux. Upon decreasing temperature below its maximum at 400 K 
A decreases. However no data have been reported below ~ 268 K. Using molecular dynamics 
computer simulations of water down to T = 220 K, we find the surprising result that the A displays 
displays a minimum at lower temperatures, Tmin = 250 K. Our simulations are for the TIP5P model 
of water which is known to reproduce experimental data well at ambient pressures 

We first perform equilibrium simulations at atmospheric pressure for temperatures from 
300 K down to 220 K in the NPT ensemble using the Berendsen thermostat and barostat with 
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a time step of 1 fs. We then run the equiUbrated configurations at different T in the NVE and 
NVT ensembles to generate the results from which we calculate A. To calculate A, we use two 
different methods: (i) the Green-Kubo relation between A and the energy current correlation func- 
tion ^, and (ii) the Miiller-Plathe algorithm ^. Our values for A obtained from both methods differ 

from each other by less than 15%. 

Figure 1(a) shows that A decreases with decreasing T, reaching a minimum at Tmin ~ 250 K, 
and increases upon further decrease in T. A clue to the possible interpretation of this surprising 
result is provided by the observation that the temperature at which A displays a minimum is the 
same temperature at which the specific heat Cp displays a maximum [see Fig. 1(b)]. 

In fact, a maximum in Cp occurs upon crossing the locus of maximum correlation length, 
the Widom line Tw{P) ^, emanating from the hypothesized liquid-liquid critical point Below 
Tw{P), water becomes locally more tetrahedral and the local structure of liquid water below 
Tw{P) more resembles ice Ih- Hence A should increase with further decrease in temperature 
below Tw.(P). 

Finally we calculate the thermal diffusivity Dt, related to p, A and Cp as 

Dt = -4r- (1) 
pCp 

Figure 1(c) show that Dt decreases very slowly upon decreasing T and rises sharply upon further 
decreasing T below the temperature of maximum Cp. 

In summary, we report a new anomaly of water, a minimum in the thermal conductivity. Our 
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findings are of special interest since they support the presence of a liquid-Uquid phase transition in 
water 
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Figure 1 : (a) Temperature dependence of the thermal conductivity A of TIP5P water, showing a 
minimum at Tmin ~ 250 K. Values of A obtained from the Miiller-Plathe algorithm are shown as 
filled circles and A obtained via the Green-Kubo relation are shown as filled squares, (b) Temper- 
ature dependence of the specific heat Cp showing a maximum at Tw ~ 250 K. (c) Temperature 
dependence of the thermal diffusivity Dt = A/ pCp, showing a sharp increase below Tw{P). 
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